Computing the Singularities of Rational Surfaces* 

O ■ S. Perez-Diaz, J.R. Sendra, C. Villarino 

. ; Dpto. de Matematicas 

O ; Universidad de Alcala 



7^ 



E-28871 Madrid, Spain 



(SI ■ {sonia.perez,rafael.sendra,carlos.villarino}@uah.es 

■<C ■ Abstract 

'•pi \ Given a rational projective paranietrization 'P(t, s, ts) of a rational projective 

2 ' surface 5 we present an algorithm such that, with the exception of a finite set 

(maybe empty) !B of projective base points of "P, decomposes the projective 
parameter plane as P^(]K) \ 5B = U^^-^6fc such that if (to : Sq : Dq) G Syt then 

CN ■ P(to,So, Do) is a point of S of multiplicity k. 

(N 

(N ■ 1 Introduction 

^ \ The study, analysis and computation of the singular locus of algebraic varieties is an 

old but still very active research topic. The interest on the study of singularities is 
motivated by multiple reasons, being one of them their applicability; for instance, in 
geometric modeling, when determining the shape and the topology of curves (either 

K> ! planar or spatial) and of surfaces, etc. In this paper, we focus on the problem of 

H \ computing the singularities, as well as their multiplicities, of rational surfaces given 

parametricaly. 

When the algebraic variety is given as a zero set of finitely many polynomials, the 
singularities and their multiplicities can be directly computed by applying elimination 
theory techniques as Grobner bases, characteristic sets, etc. However when the al- 
gebraic variety is unirational and it is given by means of a rational parametrization, 
besides the question of computing the singular locus and its multiplicity structure, 
one has the additional problem of determining the parameter values that generate the 
singular points with their corresponding multiplicities. This, for instance, can be use- 
ful when using a parametrization for plotting a curve or a surface or when utilizing 
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a parametrization for analyzing the intersection variety of two varieties being one of 
them parametrically given. Of course, one can always apply elimination techniques 
to first provide the defining implicit polynomials of the variety, second to determine 
the singularities from these polynomials, third to decompose the singular locus w.r.t. 
the multiplicities, and finally to compute the fibre (w.r.t. the parametrization) of the 
elements in the singular locus. Nevertheless, this can be inefficient because of the 
computational complexity. 

So the challenge, in the unirational case, is to derive the singularities and their 
multiplicity directly from a parametric representation avoiding the computation of the 
ideal of the variety. The case of rational curves (both planar and spatial) has been 
addressed by several authors (see [T], |1], [S], [ID])- However, the case of rational 
surfaces has not been so extensively studied. We refer the reader to [3] where the case 
of rational ruled surfaces is analyzed. 

In this paper, we present an algorithm for computing the singularities of a rational 
projective surface from an input rational projective parametrization not necessarily 
proper (i.e., birational). More precisely, the problem we deal with is stated as: 

Problem statement 

• Given a rational projective parametrization 

P(t,S,D) = (pi(t,S,D) : ■■■:p4(t,S,D)), 

of a rational projective surface S C P'^(K), where K is an algebraically closed 
field of characteristic zero. 

• Decompose F^{K) as p2(K) = ui^i&k such that if (to : Sq : Do) G &k then 
P(to,So, Do) is a point of S of multiplicity k. 

Although abusing the terminology, we will use the following definition 

Definition 1.1. The elements in ©i are called P-simple points of iS, and the elements 
in &k, with k > 1, "P-singularities of S of multiplicity k. We refer to these points as 
affine (either P-simple or P-singular) points if Do 7^ and points (either P-simple or V- 
singular) at infinity if Dq = 0. Moreover, we represent the multiplicity of (to : Sq : Do) as 
mult ((to : So : Do)) meaning 

mult((to : So : Do)) = mult(P(to,So, Do),5)) 

where mult(A, 5) denotes the multiplicity of A G P^(]K) w.r.t. S. 

The polynomials pi, . . . , p4 are assumed to be homogeneous of the same degree and 
coprime. Therefore the parametrization P(t, s, d) induces the regular map 

V : P2(K)\<B -^ S 
a ^ V{a) 



where ® = {a G P^(K) | pi(a) = ■ ■ ■ = ^4(0;) = 0}; we call the elements in 53 the base 
points of the parametrization (see Section |2]). We will be able to decompose, as above, 

03 is either zero dimensional or empty. So, we will be missing (at most) finitely 
many parameter values in P^(K). On the other hand if 53 = 0, since S is irreducible 
and V regular, then V{¥^(K)) = S (see e.g. Theorem 2, page 57, in [IS])- Therefore, if 
OS = 0, our method will determine all singularities of S. However, if ® 7^ the method 
will generate all singularities in the dense set P(P^(K) \ 03) C 5. For avoiding this 
deficiency one may consider reparametrizing normally the parametrization, however 
this not an easy task (see [9]). We do not deal with this issue in this paper. 

Our method is based on the generalization of the ideas in [5] in combination with the 
results in [6] and [7] that perform the computations without implicitizing. Intuitively 
speaking, the method works as follows; see Section [2] for further details. First we state 
a formula for computing the multiplicity of an affine point w.r.t. an affine surface 
(see Section |3]). Then, we analyze the multiplicity of the (affine) parameter values 
of the form (to : So : 1) to later study the parameter values (at infinity) of the form 
(to : 5o : 0). In order to compute mult ((to : So : 1)) we consider the four affine rational 
parametrizations (we call them V^i, ■ ■ ■ , Vx^) generated by P(t, s, o) by dehomogenizing 
w.r.t. the first, second, third and fourth component of the parametrization, respectively 
and taking = 1. Then, we apply the multiplicity formula to (to : Sq : 1) via Vx^- This 
first attempt will classify all affine parameter values with the exception of a proper 
closed set, and hence with the exception of finitely many component of dimension 
either 1 or 0. By using consecutively 7^X3,7^X2 and Vx^ we achieve the multiplicity of 
all affine parameter values not covered by Vx^^ and not being base points (see Section 
H]). Finally we deal with the parameter values at infinity with a similar strategy but 
dehomogenizing with either t = 1 or s = 1. 

The structure of the paper is as follows. In Section [2] we introduce the notation 
as well as the general assumptions that essentially imposed that S is not a plane. 
In Section [3] we state the multiplicity formula, we develop a method for computing 
a point not on the surface (this will be needed in the algorithm), starting from the 
parametric representation and without implicitizing. Moreover, we briefiy recall some 
procedures from [B] and [7]. Sections S] and deal, respectively, with the affine V- 
singularities and the "P-singularities at infinity. Section [6] summarizes all the ideas by 
deriving an algorithm. Also, a complete example is given. Sections H] and |5] contain the 
technicalities of the theoretical argumentation of our method. A reader not interested 
in that theoretical argumentation might skip these sections to go to SectionElto directly 
apply the procedure. 



2 Notation, general assumptions and strategy 

In this section we introduce the notation that will be used throughout the paper, as 
well as the general assumptions. 

IK is an algebraically closed field of characteristic zero, and P^(]K),P'^(]K) are the 
projective plane and projective space over K, respectively. Let {xi : x^ : x^ : x^) be 
the projective coordinates in P^(K) 

Let S C P'^(]K) be a rational projective surface different of the planes Xj = 0. This 
is not a loss of generality because, in that case, the surface is smooth. In addition, let 
P(t, s, d) be a rational projective parametrization of S. We consider that 7^(1,5, o) is 
expressed as 

P(t,S,t3) = (pi(t,S,t)) :p2(t,S,0) :p3(t,S,0) :p4(t,S,D)) 

where gcd(pi,p2, p3,p4) = 1 and the four polynomials are homogeneous (note that 
none of them is zero) of the same degree. We say the (to : Sq : Oo) G P^(K) is 
a (projective) base point of P(t,s,o) if pi (to, Sq, Do) = p2(to,So,Oo) = p3(to,So,t)o) = 
p4(to,So, Oo) = 0. We denote by 03 the set of (projective) base points of ■p(t,s, o). 
Since 53 is the intersection of the projective curves defined by pj(t, s, d), and since 
gcd(pi,p2, p3,p4) = 1, we get the following lemma. 

Lemma 2.1. Card(Q3) < oo. 

Our strategy will reduce the problem to the affine case to afterwards analyze the 
points at infinity. We denote by 5a;. the affine surface obtained by dehomogenization 
of S with Xi = 1; note that S is not the plane Xi = 0. Also, we denote by Vxi(t,s) 
the corresponding affine parametrization obtained from V{i,s, o). More precisely (say 
i = 4, and that 7 = deg(pj)), V{i,5, d) can be replaced by (note that p4 7^ 0) 

Pl(t.5.P) _ p2(t,5,p) p3(t,5,p) ^ 
p4(t,S,t)) ■ p4(t,S,0) ■ p4(t,S,0) ■ 

P^Pi(^^l) . t^"P2(^f,l) . ^'PsillA) . ;\ _ / Pi(^t:l) P2(^f,l) . P3(^,f^i: 



.7p,(i, £, 1) • ,7p,(i, i, 1) • o7p,(i, i, 1) • ; V P4(^ : f , 1) ' Pdl f , 1) ■ P4(i, f , 1) 
So, introducing the notation {t = - , s = -} and Pi{t, s) = pi{t, s, 1), we have that 

T> (+ \ fP^it,s) P3{t,s) PA{t,s)\ fPl{t,s) P2{t,s) Pi{t,s] 



Pl(t,s)'pi(t, S)'pi(t, S)y' ' ' \P3it,s)' P3{t,s)' P3{t,s) 

p (t s) = {M^lA ^3(^'^) P4(t,s) \ ^ (t S) - fP^^^'^") P2{t,s) P3{t,s) 



P2{t,s)' P2{t,s)' P2{t,s)J' * ' \p4{t,s)' P4{t,s)' P4{t,s) 

Observe that, since gcd(pi,p2,p3,p4) = 1, then gcd{pi,p2,P3,p4,) = 1. Further- 
more, note that Vxi {t, s) is a rational parametrization of the affine rational surface 



Sxi- Analogously, we say that (to, sq) G K^ is an (affine) base point of P(t, s, o) if 
pi(to, So, 1) = p2(to, So, 1) = p3(to, So, 1) = p4(to, So, 1) = 0. Let us denote by 53a the set 
of affine base points of V{t,s, o). Observe that Q3a can be naturally embedded in 03. 

Furthermore, we will consider that the rational functions in VxXti s) are expressed 
in reduced form. So, in the sequel, we also use the following notation 

-n (+ \ /P2(t,s) Pz{t,s) p^{t,s)\ f Pi{t,s) P2{t,s) PA{t,s) 



gi,i(t, s)' gL2(t, s)' gi,3(t, s)/' '"' ' V?3,i(t, s) ' g3,2(t, s) ' g3,3(t, s 

-p (+ ^ _ / Pi(^^^) Psjt^s) P4{t,s) \ f Pi{t,s) P2{t,s) P3{t,s) 



.92,i(t, s)' g2,2(t,s)' g2,3(t, s)y' ' V?4,i(t,s)'g4,2(t, s)' g4,3(t, s) 

where all rational functions are in reduced form. As shown in the next lemma, the 1cm 
of the denominators generates QS^. 

Lemma 2.2. 03^ = n?=i{(^o,so) e K^ | lcm(gj,i,gi,2, g^.s) = 0}. 

Proof. Let Aj := lcm(gj i, gj 2, ^i.s)- We prove that Aj = pi, and from there the proof is 
trivial. Let us assume that Ai ^ pi, similarly for the others. Clearly Ai divides pi. So, 
there exists a non-trivial factor Hi such that pi = AiHi. Therefore, by construction. 
Hi divides pk, with k = 2,3, 4. But this is a contradiction, since gcd(pi, . . . , ^4) = L n 



Furthermore, if : K" — >■ IC" is a rational affine map, we denote by deg(G) the 
degree of the map 6 (see e.g. [I3] pp.143, or [2] pp.80). In particular, deg(Pa;.) denotes 
the degree of the rational map Vxi induced by the rational parametrization VxX^-, s). 

Also, for a rational function x 'we denote by Numer(x) the numerator of x when 
expressed in reduced form. By PrimPart{„^^...^„„}(/) and Content {„^^...^t,„}(/), where 
/ G K[xi, . . . ,Xm\[vi, . . . ,Vr^, wc dcuotc the primitive part and the content w.r.t. 
{f 1, . . . , Vn} of /, respectively. For /, g polynomials depending on s we denote ReSs(/, g) 
its resultant. 

For P G P^(]K), we represent by mult(P, iS) the multiplicity of P on S; note that if 
P can be seen in the same affine space as Sx^ then mult (P, 5) = mult(P, iS^;-). 

General assumption 

We assume that for every two different polynomials pi, pj it does not exist A G K 
such that pi = Xpj. Note that if pj = Apj, then S is the plane of equation Xj — Xxj = 0, 
and the problem is trivial. In addition, note that this requirement implies that none 
of the dehomogenizations Sx^ is empty (pi 7^ and hence S is not the plane x.j = 0). 
Moreover, Sxi is not a plane parallel to any of the affine coordinate plane in K^. So 
this does not imply any loss of generality. 

Strategy 



We briefly describe here the ideas of our strategy. The precise details on how to 
execute them will come in the subsequence sections. The main steps in our strategy are 
as follows; we recall that our goal is to decompose P^(K) such that for (to : Sq : Vq) G 
P^(]K) we know whether V(tQ, Sq,Vq) is singular or simple in S, and if it is singular we 
also want to determine its multiplicity. 

1. First we analyze the parameter values of the form {to : Sq : 1) . For that, we work 
with Vx4,{t,s) and we treat the problem in K^. At this stage, we will be able 
to give an answer for K^ \ {(to, Sq) G K^ | lcm(g4 i, ^4 2, ^4,3) = 0}. Repeating the 
process (not necessarily for the whole affine plane, but for those parameters values 
that are not yet under control) with Vx^it, s), and if necessary with Vx2(t, s) and 
Vxi (t, s) we will be able to give an answer for (see Lemma 12. 2p 

4 
5^' \ n^^^O' ^0) ^ ^' I lcm(gi,i, %2, qi,3) = 0} = K^ \ ^,. 

i=l 

2. We analyze the case of (0 : 1 : 0), checking first whether (0 : 1 : 0) G 25. 

3. We analyze the case of the parameter values (1 : A : 0). First we find those A 
values generating base points. Afterwards, we study (under a suitable dehomog- 
enization) the rest of the points. 

3 The multiplicity formula 

In this section we state a formula for computing the multiplicity of a point in K'^ 
w.r.t. an affine rational surface in K?, when a rational parametrization (not necessarily 
proper) is provided. As a consequence, we give a criterion for an affine point to be 
regular on the affine rational surface. In order to derive an algorithmic version of these 
results, we will recall some procedures in [0],[7] and [5], and we will present a method 
for determining a point out of the surface without knowing the implicit equation. 

For that purpose, throughout this section, 2 C K^ is a rational affine surface and 

0(t ) = (MA ^2(t,g) Nsjt.s) 
^^''^ \D^{t,syD2{t,sy Dsit,s) 

a rational parametrization (in reduced form) of Z; we assume w.l.o.g. that Z is not a 
plane parallel to the coordinate planes of K^. 

For any A = (a, b, c) point of K^, let /(x, y, z) be the defining polynomial of Z and 
F(x, y, z, w) its homogenization. We consider the polynomial g{x, y, z) = f{x + a,y + 
b,z + c), as well as G{x,y, z,w) = F{x + aw,y + bw, z + cw,w). It is clear that 

mult(A,Z) = deg(G)-deg^(G). 



On the other hand, note that 

Ni{t,s) N2{t,s) Ns{t,s) 

parametrizes the projective surface defined by G. Therefore, since we have assumed 
that Z is not a plane parallel to the coordinate planes, then Ni/Di — a ^ and hence 

^^ N2it,s)-bD2it,s) Di{t,s) Ns{t,s)-cD3{t,s) Di{t,s) Di{t,s) 



Ni{t,s) -aDi{t,s) D2{t,s)' Ni{t,s) -aDi{t,s) D3{t,s)' Ni{t,s) -aDi{t,s] 

parametrizes the affine surface defined by G{l,y, z,w); note that, since G is homoge- 
neous, deg^(G) = deg^{G{l,y, z,w)). Let us introduce the following notation 

^ ^ r N2{t,s)-bD2{t,s) Di{t,s) Ns{t,s)-cD3{t,s) Di{t,s) 

'''^ ^^'^^ \Ni{t,s)-aDi{t,s)' D2{t,sy Ni{t,s)-aDi{t,s)' D3{t,s) 

and let $2,3(^) : K^ -^' K^ be the induced map. Moreover, if for i = 1,2, xf(t, s) 
denotes the i-component of $2,3(^)(i, s), let g^ ' be the polynomial 

gf '^(t, s, hi, /12) = Numer(xf (t, s) - xt{hi, 112)), i = 1,2 

where hi, /i2 are new variables, and let K{t, s, hi, /i2) = gcd(g]^ ' , 02 ' ) where the gcd 
in computed in K[hi, h2][t, s]. Then we introduce the polynomial 

nQ'Af. , > h \- ! K{t,s,hi,h2) if deg|t^,}(K) >0 
[t,s,ni,n2)-<^^ ifdeg^,,j(K)=0 

Remark 3.1. We observe that if 0^'"^ = 1 then the determinant of the Jacobian of 
$2 3(^)(^, s) is not identically zero (see the preliminary paragraphs to Theorem 1 in 

In the following theorem and corollaries we assume that: 

1. none of the projective curves defined by each of the non-constant polynomials in 

{A^i, N2, N3, Di, D2, D3} passes through (0:1:0) 

2. for each A = (a, b, c) G K^ (similarly Aq) none of the projective curves defined 
by the each of the non-constant polynomials in {A'^2 ~ bD2, N^ — cD^, Ni — aDi} 
passes through (0:1:0). 

Note that, if necessary, one can always perform a suitable polynomial linear change of 
parameters. 

Theorem 3.2. (The general formula) It holds that 
1. mult(A^) = deg(Z) iffQ^'^ ^ 1. 

7 



2. mult{A,Z) < deg(Z) iff Q^'^ = 1. Furthermore, if Q^'"^ = 1 then 

deg(<l>23(A)) 



deg(Z) -mn\t{A,Z) 



deg(Q) 



Proof (1) By Theorem 4 in [7], deg^(G) = iff fl®'^ ^^ 1. Hence, mult(A, Z) = deg(Z) 

iff flfi.A _^ 1^ 

(2) All liypotliesis of Tlieorem 6 in [S] are satisfied. Tlierefore, deg^(G) = "fe (q*) ' 
and tlie result follows by taking into account that deg(Q) = deg(Q*). n 

Remark 3.3. Note that: 

1. If 02,A _^ i^ by Theorem 0(1), then AeZ. 

2. For A ^ Z, deg($2,3(^)) is invariant, indeed it is deg(Z)deg(Q). 

3. Let A = (a, b, c) and assume that a 7^ 0; similarly for b and c. We embed 
A in P^(K) as (a : fe : c : 1). Let 2 be the projective closure of Z, and 
let Z^-^^ the dehomogenization of Z w.r.t. Xi = 1. Let Q^,'^ be the corre- 
sponding parametrization of Z^-^ generated by Q. Then deg(Z) = deg^Z^J, 
mult((a, 6, c),Z) = um\t{{b/a,c/a,l/a), Z^J and deg(Q) = deg(Qa;J. There- 
fore, deg(<l>2,3((a,^c))) = deg(<l>2,3(Vo>c/a, 1/a))). 

Corollary 3.4. Z is a plane if and only if there exists a non-empty dense subset Q of 
Z such that for each A eD., q^'^ 7^ 1. 

Proof. Since Z is irreducible, the corollary follows from Theorem 13.21 (1). n 

Corollary 3.5. Let Z not be a plane. Z is a cone of vertex A if and only if g^'"^ 7^ 1. 

Proof. Let Q{t,s) = t{gi{s) , g2{s) , g^ls)) be a rational parametrization of Z where 
we assume w.l.o.g. that A is the origin O. The defining polynomial f{x,y,z) of Z 
is a form of degree d > 1. Therefore mult(0,2) = d = deg(^) and by Theorem 
K2\ gS^o _^ i_ Conversely, if q^'^ ^ 1 then deg^(G) = (see Theorem 4 in [7]). 
Thus, G{x,y,z,w) = g{x,y,z) is an irreducible form of degree d > 1. Let us see that 
{g{x,y, z) = 0} n {x = 1} defines a rational plane curve V of degree d. Since g is a. 
form and Z is not a plane, D is a curve of degree d. Moreover, since ^ is a surface, 
^2,3{A){t,s) is not constant and parametrizes the surface defined by g{l,y,z). Then 
substituting either t or s by a suitable constant (say t by to), T^{s) := (1, $2,3(A)(to, s)) 
parametrizes V. 

Now, Q{t, s) = A + tTZ{s) defines a cone of vertex A contained in Z. Therefore, 
since Z is irreducible, it holds that Z is the previous cone. n 

Corollary 3.6. Let Z not be a plane. There exists at most one A G K^ such that 



Proof. Let us suppose that exist two points verifying the statement. We assume 
w.l.o.g., that one of them is the origin O and the other one is A = {a,b,c) ^ O. 
By Corollary I3.5[ Z is a cone parametrized as Q(t, s) = t{gi{s) , g2{s) , g3{s)) . Note that 
f{x, y, z) is homogeneous of degree d > 1. Moreover g{x, y, z) = f{x + a,y + b,z + c) 
can be expressed as g{x, y, z) = f{x, y, z) + g*{x, y, z) where deg((7*) < d. By Theorem 
I3.2[ mult(y4, Z) = deg{Z). Thus f{x + a,y + b,z + c) = f{x, y, z). This implies that 
Z is invariant under the translation of the vector A. Let us see that 2 is a plane 
which is a contradiction. Indeed, let Sq G IK be such that v := (5'i(so)!5'2('So),fi'3(so)) is 
not parallel to A; observe that sq exists because {gi{s) , g2{s) , g2,{s)) is not a line. We 
consider the plane LI given by the parametrization p{s,t) := sA + vt. Because of the 
invariance of Z, under translation of vector A, the family of lines p{n, t) are included 
in n n ^, for n G N U {0}. Since all these lines are different and Z is an irreducible 
algebraic set one concludes that Z = 11. n 

Corollary 3.7. (The multiplicity formula) Let Aq eK^\Z and let A G K^. Then 

1. ifQQ'^ = 1, then 

raum,Z) = deg('l'.3(^))-des(t.3(^))^ 

deg(Q) 

2. ifQ^'^^ 1, then 

^ ' deg(Q) 

Proof. We prove (1); similarly for (2). By Theorem 13. 2| one has that {deg{Z) — 
mult(A, Z))deg(Q) = deg($23(^)), and deg(Z)deg(Q) = deg($23(Ao)). From here the 
proof is obvious. n 

Remark 3.8. Note that 

1. if there exists Aq such that g^.A) -^ \ {}.e., Z is a. cone), Theorem 13.21 and 
Corollary 13.71 provide the degree of the surface. 

2. From Corollary 13. 7[ one deduces that deg($23(^)) is invariant for all the A G K'^ 
having the same multiplicity w.r.t. Z. 

The next corollary is a direct consequence of Corollary 13.71 

Corollary 3.9. (Criterion for simple points) Let Aq E K^ \ Z and let A G K^. If Z is 
not a plane, the following statements are equivalent 

1. A is a simple point of Z. 

2. gQ'^ = 1 and deg($2,3(^o)) - deg($23(^)) = deg(Q). 

9 



Proof. By Theorem [321 (deg(Z) - l)deg(Q) = deg(<l'2,3(^)), and deg(Z)deg(Q) = 
deg($2,3(^o))- From here the proof is obvious. n 

We observe that if we know how to compute deg(Q), deg($23(^)) for any given 
A G K^, and if we know how to compute a point out of the surface (recall that we do 
not have the implicit equation of Z), Corollary 13.71 provides a method for computing 
the multiplicity of any point in K.^, and Corollary 13.91 a method to check whether it is 
simple on the surface. We note that the deg(Q) is the index of improperness of Q(t, s); 
if Q{t, s) is proper then this index is 1. Therefore, once the parametrization is given, 
deg(Q) is fixed. However, deg($23(^)) will vary depending on A. Both quantities 
can be derived by applying elimination theory techniques as Grobner basis. Indeed, 
they can be computed by means or resultants as shown in [6] without determining the 
implicit equation of Q. 

In the following we recall (as a recipe) how to compute deg(Q) and deg($23(^)); 
for further details we refer to [B], [7] or [S]. In addition, we deduce a method for 
determining a point out of the surface. 

Method 1: Computation of deg(Q) 

[Step 0] Check the global hypothesis above: 

[Step 0.1] If any of the projective curves defined by the each of the non-constant 
polynomials in {Ni, N2, N3, Di, D2, D-s} passes through (0:1: 0), apply a suitable 
(polynomial) linear change of parameters. 

[Step 0.2] If the determinant of the Jacobian of (^, ^) is identically zero, apply 
a suitable linear change of coordinates in K^; namely, exchange suitably the affine 
coordinates in K^. 

[Step 1] For ^ = 1, 2, 3, compute G,{t, s, h, h^) = Numer (fM _ §|iM) . 
[Step 2] Determine /2(t, /ii, /i2, X) = ReSs(Gi, G2 + ^^3) where X is a new variable. 
[Step 3] Compute S(t,hi,h2) = PrimPart{/jj^/i2}(Contentx(-R)))- 
[Step4]deg(Q)=deg,(^). 

Computation of deg($23(^)) 

We observe that by Theorem 13. 2[ ii Aq ^ Z then g^'"^" = i^ Therefore, by Corollary 
13. 7[ we only need to compute deg(<l>23(^)) for those A such that q^'^ = 1; in particular 
when A = Ao- Thus, in the following we assume that A = (a, b, c) is such that 
qQ,a _ ]^ Moreover, we will use the following technical lemma that will simplify the 
computations. 

In addition, since g^'^ = 1, by Remark 13. 1[ the determinant of the Jacobian of 
^2,3{A)it,s) does not vanish. Therefore, $2,3(^)(K^) is dense in K^. So, we can 
compute the degree by taking a generic element (Xi,X2) G K^ as it is done in [8]. 
More precisely, we have the following method. 
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Method 2: Computation of deg($23(^)) 

[Step 0] Check the global hypothesis above: if any of the projective curves de- 
fined by the each of the non-constant polynomials in {A^^i — aDi,N2 — bD2,N3 — 
cD^, Ni, N2, N3, Di, D2, Ds} passes through (0:1: 0), apply a suitable (polynomial) 
linear change of parameters. 

[Step 1] We take the components xfit, s) of $2,3(A)(t, s) := {xiit, s),x^(^, s)) 
[Step 2] For i = 1, 2, let Hi{t, s, X^) = Numer(xH*. s) - X») G K[t, s, X^]. 
[Step 3] i?(t, Xi,X2) = ReSs{Hi.H2) E K[t,Xi,X2] 
[Step 4] deg($2,3(^)) = deg,(PrimPart{x„x,}(^)). 

Computation of a point Aq out of Z 

For our reasoning we need to know the partial degree, w.r.t. one of the variables, 
of the defining polynomial of Z. Say that m is the partial degree w.r.t. the variable x 
(below we show how to compute m). This means that for almost all affine lines C of 
the type {y = X, z = fi} (recall that Z is not a plane parallel to the coordinate planes) 
it holds that Card(£ (1 Z) = m. Then, the idea is as follows. We take values for (A, /i) 
till the number of different points on Z generated by Q{t, s) is m. Note that for a fixed 
(A,yu), these points are 






D2{to,So)X-N2{to,So)=0, 

Dsito, So)/" - iV3(to, So) = 0, 
lcm{Di,D2,Ds){to,so)y^O 



Once we have found a suitable (A,/i), every point (a, X, fi) ^ >V(A,/i) is not on Z. We 
finish this section showing how to compute m (see details in Theorem 6 in [8]). 

Method 3: Computation of the partial degree w.r.t. x 

[Step 1] Apply Method 1 to compute deg(Q) 

[Step 2] Let, for i = 2,3, Gi(t, s, hi, /i2) as in Step 1 of Method 1. 

[Step 3] Return - — —degt{FnmPaTt{huh2}(^^^s{G2{t,s,hi,h2),G3{t,s,hi,h2))). 

Remark 3.10. Note that the polynomials Gi are obtained in Step 1 of Method 1, 
and therefore it might happen that Step of Method 1 was required. In that case, 
we would have performed a linear change in the parameters {t,s}, and/or an affine 
linear change of coordinates {x, y, z} consisting in a permutation of variables. The 
first situation does not affect to the partial degree of the polynomial. However, the 
second can. Nevertheless, if this is the case, we only need to work with the new 
variable (the one exchange with x) and the corresponding lines perpendicular to its 
corresponding coordinate plane. 

Method 4: Computation of Aq eK^\Z 
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[Step 1] Apply Method 3 to compute the partial degree m of the defining polynomial 

of Z w.r.t. X. 

[Step 2] Give values to (A,/i) G K^ till Card(W(A,/i)) = m, then take Aq := (a, A,/i) G 

K3\W(A,/x). 

Computation of um\t{A, Z) 

We finish the section, putting together all the previous ideas for computing the 
multiplicity of A G K'^ w.r.t. to the rational affine surface Z, parametrized by Q{t, s). 

Method 5: Computation of mult (A, Z) 

[Step 1] Apply Method 4 to find a point Aq ^ Z. 
[Step 2] Compute q^'^. 
[Step 3] If gS.A -^ I then 

[Step 3.1] Apply Method 2 to compute rii := deg($2,3(^o))- 

[Step 3.2] Apply Method 1 to compute n2 := deg(Q). 

[Step 3.3] Return ^ 
[Step 4] If 02,A = I then 

[Step 4.1] Apply Method 2 to compute rui := deg($2,3(^)) and rii := deg($2,3(^o))- 

[Step 4.2] Apply Method 1 to compute n2 := deg(Q). 

[Step 4.3] Return ^^i^ 

4 Computing the affine ^-singularities 

In this section we see K^ embedded in P^(K) by means of the natural map 

j:K2^p2(K),(to,So)^(to:so:l); 

in this sense, as already commented in Section [T] we will be determining the affine 
P-singularities of S. 

For this purpose, let Aj := {(to, Sq) \ lcm(gj i, qi2, ?i,3)(to, Sq) = 0} and 55^ be the set 
of base points of Vxi(t,s). Note that j(23a) C 55. The basic idea consists in applying 
Method 5 to a generic point on S. For this purpose, we proceed as follows. 

First Level. We decompose Ai := K^ \ A4 as 

such that if (to, Sq) G ^l then Vx^ito, Sq) is a point of S.^^ of multiplicity k. 
Second Level. If A4 \ Q3a 7^ we decompose A2 := A4 \ A3 as 

A2 := ^iUtu 
such that if (to, Sq) G ^1 then ^^3(^0, ■^o) is a point of S^^^ of multiplicity k. 
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Third level. If A2 \ 23a 7^ we decompose A3 := (A4 fl A3) \ A2 as 

A3 = uiudl 

such that if (to, sq) G ^l then "P^^l^o, sq) is a point of Sx2 of multipUcity k. 
Fourth Level. If A3 \ ®„ 7^ we decompose A4 := (A4 fl A3 n A2) \ Ai as 

A4 = uiudt 

such that if (to, Sq) G Jt then Vx^it^, Sq) is a point of S^^ of multiplicity k. 

Note that at this point, A4 \ ^„ = 0. Moreover }{V^U^'^{) C &k (see Section [1]). 

First level 

The strategy for this level is as follows. We determine a closed set A* of K^ such 
that for every (to, Sq) G Ai\ A* then Vx^^ito, sq) is simple on Sx^] note that Ai\ A* C 5i- 
Next we decompose Ai fl A* as 

Ai n A* = U^f^^^A{k) 

such that if (to, sq) G A(A;) then mult(P(to, sq, 1),S) = 4- Note that A(A;) C d}^- 

First level (Part I): computation of A* 

In order to compute A* we will determine some closed sets {A*}j=o,..,4 of IK^ such 
that A* = Uf^^A*. For that purpose, we apply Method 5 in Section [3] to a generic 
point of Px4(Ai) C Sx^] namely Vx^(t*,s*) where t*,s* are treated as new variables. 
We assume that we have already computed a point Aq in ]K'^\5^^ (see Step 1 in Method 
5 or see Method 4 in Section [3]) as well as deg(<I>2,3(^o)) ^"^^ ^^sO^x^)- For simplifying 
the notation, throughout this section we will denote the generic point 7^X4 (t*,s*) by 
A*. 

To perform Step 2 in Method 5, we consider ^2,3{A*){t,s) as well as its rational 
function components; namely 

A*U„^* „*N _ 54,2(t*,S*)p2(t,s)-p2(i*,S*)g4,2(t,s) qi,l{t, s) q^,l{t* , S*) 






g4,i(t*,s*)pi(t,s) -pi(t*,s*)g4,i(t,s) g4,2(t, s) g4,2(t*,s*)' 
g4,3(t*, s*)p3(t, s) - p3(t*, s*)g4,3(^, s) g4,i(t, s) qA,l{t* , s*) 



g4,i(t*, s*)pi{t, s) - pi{t*, s*)g4,i(t, s) q4,3{t, s) q4,3{t*, s*) 

Note that, since t, s, t*, s* are independent variables, and since we have excluded planes 
parallel to the coordinate planes (see general assumptions in Section [2]), the above 
rational functions are well-defined. Moreover, for every particular value (to, sq) G Ai 
of (t*,s*) the specialization of the rational functions are also well-defined. Similarly, 
we take the polynomials 

gf^--^* = Numer(xf (t, s, t*, s*) - xf {h, h^, t\ s*)) 
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as well as 

K{t,s,hi,h2,t*,s*) =gcd{Q^"^' ,92^*' ), 

and g'^^4'-4 . That is, we perform Step 3 for the generic element A*. 

If g^^4"^* ^ 1, we do not need to continue since, by Corollary 13 ■4[ 5 is a plane and 
the problem is trivial. Alternatively, one might avoid this case by trivially checking 
first whether 5 is a plane. Therefore, we assume that g^'^i'-^* = 1. Thus, we are already 
in Step 4 of Method 5. Nevertheless, at this stage, we know that generically g^^4'-4* = i 
but for certain (t*, s*)-values the gcd may increase the degree. Note that, because of 
Corollary 13. 51 and 13. 6[ this can only happen if iS is a cone and only for the (t*, s*) values 
that generate its vertex. These values, if they exist, will be included in the closed set 
Aq. For determining Aq we consider the following direct generalized version of Lemma 
3 in [H]: 

Lemma 4.1. (Lemma 3 in [11]) Let f,g E K[t* , s*][Y][t] \ {0}, where Y is a finite set 
of variables. Let f = f ■ gcd{f,g), g = g ■ gcd{f,g). Let {to, sq) G K^ be such that not 
both leading coefficients of f and g w.r.t. t vanish at (to,So). IfKest{f,g) does not 
vanish at (to,So), then gcd{f, g){t,Y ,to, sq) = gcd{{f{t,Y ,to, sq), g{t,Y ,to, sq)). 

In our case, we consider q^ ""■*' G K[t*, s*][/ii, /i2, s][t]; note that both polynomials 
are not identically zero because we already know that S^^ is not a plane. Let Tj be 
the leading coefficient of g^ ''''' w.r.t. t, and let 

r(r,s*,/ii,/i2,s) := Resi(gi"'" ,52""" )• 

In addition, since we have assumed that g^^4'-^* = 1^ we know that K G K[/ii, /i2, t*, s*]. 
Let Zk be the zero set of at least one non-zero coefficient, w.r.t. {hi,h2}, of the 
homogeneous form of maximum degree of K. 

We define Aq as the zero set of all coefficients of Ti, T2 w.r.t. {hi, /12, s} union the 
zero set of all coefficients of r w.r.t. {hi, /12, s} union Z^- 

Now, we proceed with Step 4 of Method 5. We have assumed that ni and n2 (in 
Step 4 of Method 5) have been already computed. So, it only remains to analyze the 
determination of mi := deg($2,3(-^*))- Therefore, we apply Method 2 to $2,3('^*)(^7 s). 

We assume that none of the projective curves defined by the non-constant polyno- 
mials in {pi, P2, P3, Qa,!, (l4,2, (l4:,3} passes through (0:1: 0). If this is not the case, we 
perform a suitable polynomial linear change in the parameters {t, s}. Note that, in 
this situation, $2,3("^*)(^, s) satisfies the conditions in Step of Method 2, seeing the 
projective curves in P^(K(t*,s*)) where K(t*,s*) is the algebraic closure of ]K(t*,s*). 
However, it might happen for some particular values of {t*, s*} then condition fails. In 
order to control this, we introduce the following set Z^d- We take the homogenization 
(in the variables {t, s}) of the numerators and denominators of xf* (ty s, t* , s*), and we 
substitute them in (0 : 1 : 0). Observe that, as remarked above, the resulting polyno- 
mials are not identically zero. Now, Zoo is the union of the zero sets in K^ of these 
polynomials. 
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In Step 1 of Method 2, we take Xi *' ■, ^ = 1)2, and in Step 2 of Method 2, we 
compute 

i7,(t,s,X„r,s*) = Numer(xr^^'-^* -X,) G K[r, s*,Xi,X2, t,][s]. 

For z = 1,2, let Mj(t, Xi,X2, t*, s*) be the leading coefficient oi Hi w.r.t. s. Then, we 
define AJ as the zero set of all coefficients of Mi w.r.t. {t,Xi,X2} union the zero set 
of all coefficients of M2 w.r.t. {t, Xi,X2} union Z^ . 

In Step 3 of Method 2, the resultant polynomial R is computed. We observe that since 
qTx^A* = 1^ 7? is not identically zero. We see i? as a polynomial in K[t*, s*][t, Xi,X2], 
and hence we denote it by i?(t,Xi,X2, t*, s*). Let iy(Xi,X2, t*, s*) be the leading 
coefficient of R w.r.t. t. Then, we define Ag as the zero set of all coefficients of W 
w.r.t. {Xi,X2}. 

In Step 4 of Method 2, first we express i? as a polynomial in {Xi,X2} as 

where we collect the non-zero coefficients of R w.r.t. {Xi, X2}. We want to control the 
behavior of the primitive part under specializations, which essentially means to control 
the content. More precisely, let 

a{t,t*,s*) = gcd({aij | (z,j) G J}) = Content{Xi,x2}(^), 
and let _ 

'^^ ' a{t,t*,s*) 

Let N{t*, s*) be the leading coefficient of a w.r.t. t. We analyze (under specializations) 
the gcd of {ajj | (z, j) G J}. We distinguish several cases depending on the cardinality 
of J; we observe that Card(J) 7^ 1 since deg($2,3('Pa:4(^0) ■^o))) > 0. 

[Case 1] Let Card(J) = 2; say J = {{io,Jo),{h,ji)}- We apply Lemma WA] (i.e.. the 
adaptation of Lemma 3 in [11]) to ai^^j^jTii-^j-^, seen as polynomials in ]K[t*, s*][t]. Let 
Loit*, s*) be the leading coefficient oiTEi^j^ w.r.t. t, Li{t*, s*) be the leading coefficient 
of Oj^j^ w.r.t. t, and let S(t*,s*) = ReSi(ajojo, aj^jj. Then, we define A3 as the zero 
set of {Lq, Li} union the zero set of S, and A4 as the zero set of A^ (see above). 

[Case 2] Let Card(J) > 2; say J = {{ik,jk)}k=i,...,e, with £ > 2. We apply Lemma 9 in 
[6]. For convenience of the reader we recall here the part of that lemma that we will 
use. 

Lemma 4.2. (Lemma 9 in 0) Let /, G K[r,s*][t] \ {0}, /, = /, ■ gcd(/i, . . . , /^), 
i = 1, ... ,771. Let (to, So) £ K^ be such that the leading coefficient of fi w.r.t. t does not 
vanish at (to, Sq). //ReSi(/i, /2 + EII3 ^i-2fi){to, So) 7^ 0, where Wj, j = 1, . . . , m-2, 
are new variables, then gcd{fi, ..., frn) {to, Sq) = gcd(/i(to, Sq, t), . . . , /^(to, So, t)). 
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Thus, we apply the lemma to {ai,^jf,}k=i,...,e seen as polynomials in K[t*, s*][t]. Let 
L(t*,s*) be the leading coefficient of Oj^j^ w.r.t. t, and let 

/ ^ 

S{t*, s*, Wi,..., Wt-2) = Rest ( ai^,h,ai2,j2 + ^ Wk-2ai^,]u 

\ k=3 

We define A3 as the zero set of all coefficients of S w.r.t. {Wi, . . . , W£-2}, and A^ as 
the zero set of L union the zero set of A^. 

Note that, since S is irreducible, and Vx^it, s) is a generic element of Sx^, we have 
the following lemma. 

Lemma 4.3. Let (tc-so) ^ ^1 be such that 

R(t,Xi,X2,t*, s*] 



deg($2,3(^x4(io,so))) = degt 



a(t,t*,s*) 

It holds that Vx^^ito, Sq) is a simple point of Sx^- 

We finish this subsection with the following theorem. 
Theorem 4.4. V(to, sq) G Ai \ A*, Vx^itQ, sq) is a simple point of Sx^- 

Proof. Let (to, s^) G Ai \ A*; throughout the proof, we denote Vx4,{to, sq) by Aq. Since 
(to, -So) G Ai, then (to, so) ^ A4, and hence Aq is well defined and it is a point on Sx^^- 
Moreover, Xi "{ti s) are also well-defined. On the other hand, since (to, So) ^ Aq, then 
r(to. So, hi, /12, s) does not vanish and at least one the polynomials Ti(to, Sq, hi, /12, s), 
'^2{to, So, hi, /i2, s), does not vanish. Then, by Lemma [4. 1|, 

K{t, s, hi, h2, to. So) = gcd(0r^'*-^«, gj^--^"). 

Furthermore, since (to, So) ^ Aq, then (to, Sq) ^ Z^, and hence 

deg{t,s}(^(^, s, hi, h2, t*, s*)) = deg^t, ^^{K{t, s, hi, /is, to, So)). 

Therefore, 

Q^-^'-^o^t, S, hi, h2) = S^^*'-^*(t, S, hi, h2, So, to) = 1. 

Note that since (to,so) ^ A^, then (to, so) ^ Z^ and hence the conditions in Step 0, 
Method 2, are satisfied. Moreover, neither Mi(t,Xi,X2, to, so) nor M2(t, Xi,X2, to, so) 
vanish. Similarly, since (to,so) ^ Ag, iy(Xi,X2, to, so) does not vanish. 
If we are in case 1, since (to, Sq) A3 we get that Lo(to,So) 7^ or 
Li{to,So) 7^ 0, and S'(to, So) 7^ 0. Thus, by Lemma Wl\ we get that a(t, to,So) = 
gcd(aj(,jo(t,to. So), ajiji(t, to. So)). Moreover, by well-know properties of resultants, 
we get that (up to multiplication by a non-zero constant), -R(t, Xi,X2, to, So) = 
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ReSs(ifi(t,s,Xi,to,so),-ff2(^,s,X2,to,'So)). Furthermore, since W{Xi,X2,to, so) ^ 
(see above), 

deg,(i?(t,Xi,X2,to,So)) = deg,(i?(t,Xi,X2,r,s*)). 

On the other hand, (to) -So) ^ ^4 imphes that deg^(a(t, t*, s*)) = deg^(a(t, to, -So))- 
Summarizing, 

deg($2,3(^S)) = degi(PrimPart{Xi,X2}(i?(t,Xi,X2,to,So))) = 

deg,(PrimPart{x„x,}(i?(t,Xi,X2,r,s*))) = deg(<l>2,3M'^)). 

Therefore, by Lemma [4. 3[ ^q is simple. 

If we are in case 2, since (to, Sq) ^ A3 we get that 5(^0, Sq, W\^ . . . , W^_2) 7^ 0. Since 
(to, ■So) ^ A4 we know that L(to, Sq) 7^ and A^(to, Sq) 7^ 0. Thus, by Lemma H^ we 
get that a(t, to,So) = gcd({aij(t, to, so) | (hj) G J})- From here the proof follows as in 
the case 1. D 

First level (Part II): decomposition of A* \ A4 

We decompose A* as union of irreducible closed sets; note that they are of dimension 
less or equal 1. Let C be an irreducible curve in A*. If C C A4, there is nothing to 
do. If not, we compute the intersection of A4 and C (note that A4 is empty or a plane 
curve). This intersection would be either empty or finitely many points. For an open 
subset of C, the degree of the corresponding map $2,3 would be invariant, and hence 
all points in the open subset would generate points on S with the same multiplicity. 
The complementary of this open subset is now either empty or a finite set of points. 
So, if it is not empty, we apply the formula to each of the finitely many points in the 
closed set as well as for those points in the zero- dimensional components of A*. 

In order to compute the open subset of C, we do an analogous reasoning as in the 
previous subsection. 

[Rational case] If C is rational, we compute a proper normal rational parametrization 
Q(A) of C (see [H]). Then, we apply Method 5 to n{X) := V^^{Q{X)y, say that H(A) 
is expressed as: 

-^M) _ f'fiW f2{x) ¥>3{xy 



,0l(A)' 02(A)' 03(A). 

where gcd(c/?j,0j) = l,i = 1,2,3. Note that Step 1 as well ni,n2 (in Steps 3, 4) were 
already computed in Level I (first part). In Step 2, we have to compute g'^^4'^(-'*). For 
that we distinguish two cases: 

1. if C ^ Ag (see first part of the proof of Theorem 14. 4 p then, if (to, sq) := Q(Ao) ^ 
(Cn A;;) \ A4, it holds that gP-iM^o) = 1. por the others, the finitely many 
(maybe empty) points in (C fl Aq) \ A4, one applies directly the whole Method 
5. Observe the connection with cones; see Corollaries 13.51 and 13.61 
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2. If C C Aq, we repeat the reasoning done for the computation of A* in Part 
I of Level 1st. That is, we compute g^^4.'^(-'*) generically (i.e., treating A as a 
transcendental element). Note that this, essentially, means computing a gcd in 
unique factorization domain ]K[X,hi,h2,s,t] that (see e.g. Section 4.1. in ^^) 
can be reduced to the computation in the Euclidean domain K(A, /ii, /12, s)[t]. 
For an open subset A of C, q^='4''^W = 1^ and we can go ahead through Step 4. 
For the complementary closed set (that is empty or finite) we execute the whole 
Method 5. 

Therefore, after performing the above considerations, we can assume that g^^4'^(-^) = 1. 
Thus, we pass to Step 4, and hence it only remains to apply Method 2 to compute 
deg($2,3(^(-^))! where A belongs to a non-empty open subset of K; namely those A 
such that Q(A) G C \ Aq if we come from case 1 (above) or Q(A) G C \ A if we come 
from case 2 (above). We observe that all computations can be carried out: we have 
to compute resultants in the unique factorization domain K[A, t, Xi,X2][s] and gcds in 
the Euclidean domain K(A)[t]. 

[Positive genus case] If C is not rational, we work over the field of rational functions ]K(C) 
of the curve (see [12]). Let f{t,s) be the defining polynomial of C, then K(C) is the 
quotient field of K[t, s]/(/). Then, we apply Method 5 to V(t, s), where t, s G K{C) are 
representatives of the equivalent classes oft, s respectively, i.e., t — t, and s — s belong to 
the ideal (/). We recall that the arithmetic in the field K(C) can be executed by using 
the defining polynomial of C. We observe that all computations can be carried out: 
we have to compute gcds in K(C)[/ii, /i2, s, t] (which can be performed in the Euclidean 
domain K.{C){hi, h2)[t]), resultants in the unique factorization domain ]K(C)[t, Xi,X2][s] 
and gcds in the Euclidean domain K(C) [t] . 

For each 1-dimensional component C of A* we will get an open subset where all 
points (i.e., parameter values) behave the same; that is all have the same multiplicity. 
So each of these open subsets will be part of ^l for some k. The complementary of 
these open sets are either empty or zero- dimensional. So we will have, in the worst 
case, a set of finitely many parameter values to be classified. For each of them we 
apply Method 5, and we determine their multiplicity. Finally, they are included in the 
corresponding ^J,. 

Second, third and fourth levels 

Let A4\^(j 7^ 0. We want to decompose A2 (i.e., A4\A3). We observe that A2 would 
be either empty or 1-dimensional; since Aj are either empty or plane curves. Clearly, 
the interesting case is when dim(A2) = 1. Then, for each irreducible component of A2 
we proceed as in the first level (part II). Finally, note that the same argument and 
strategy is valid for the third and the fourth levels. 



5 Computing the ^-singularities at infinity 

In this section, we show how to proceed with the steps 2 and 3 of our strategy (see 
Section |2]). So, first we analyze whether A = (0 : 1 : 0) is a P-singularity. For this 
purpose, we check whether AG®. If A ^ 03, then at least one of the polynomials pi 
does not vanish on A (say w.l.o.g. P4). Then, we replace P(t, s, o) by 



Pia,i,f).P2(|,i,f).P3(i,i,f; 



: 1 



MIA,1) P4(^l,f) P4(^l , 

So, introducing the notation {i = ^,v = ^} and pi(t, v) = pj(t, l,v), we get 

'piit.v) P2ilv) psilv) 



V.,{t,v) 



pi{t,vy Pi{t,vy Pi{t,v] 



that parametrizes Sx^- Similarly, if necessary, we introduce VxX^iv) with i = 1,2,3. 
Now, we apply Method 5 to compute 

mult((0 : 1 : 0)) = mult(P,,(0, 0),5,J = mult(P(0, 1, 0),5). 

Now, it only remains to analyze the points in £^ = {(1 : Aq : 0) | Aq € K}. For that, 
first we determine those points in S that are base points, namely 

£* = {(1 : Ao : 0) I gcd(pi(l. A, 0), P2(l, A, 0), P3(l, A, 0), P4(l, A, 0))(Ao) = 0} 

There exists i such that pj(l. A, 0) is not identically zero, since otherwise would divide 
gcd(pi, . . . , P4), which is a contradiction. Let us assume w.l.o.g. that p4(l. A, 0) is not 
identically zero. We then introduce the finite set 

r* = {(l:Ao:0)|p4(l,Ao,0) = 0}\r, 

and we proceed to compute the multiplicity of each (1 : Aq : 0) G S**. For that, we 
observe that there exists j 7^ 4 such that pj(l, Aq, 0) 7^ 0, and we apply the multiplicity 
formula using the dehomogenization of V{i,5, o) w.r.t. the j-component. 

To analyze the open subset £ \ S**, we replace P(t, s, o) by 

' Pl(l,f,^) P2(l,f.f) P3(l,f,f) 
,P4(l,f,^)-p4(l,f,l)'p4(l,f,?)- 

So, introducing the notation {s = j,v = ^} and pi{s,v) = pi{l,s,v), we get 

'piis,v) P2{s,v) Pz{s,vy 



VxAs^v) 



P4[S,V) Pi[S,V) P4[S,V) 



that parametrizes Sx^^- Similarly, if necessary, we introduce VxX^i"^) with i = 1,2,3. 
Now, one has to proceed as in Section |U Level 1 (Part II, case rational) with the 
rational curve Q(A) = (A,0). 
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6 Algorithm and Example 

In this section we summarize all the previous ideas to derive an algorithm that we 
illustrate with a complete example. For this purpose, let S C P^(K) be a projective 
surface, and P(t, s, o) a parametrization of S expressed as 

P(t,S,0) = (pi(t,S,t)) :p2(t,s,0) :p3(t,S,D) : Pi{i,s,V>)) 

where pi G K[t, s, t)] are homogeneous polynomials of the same degree, and 
gcd(pi,p2,p3,p4) = 1. Let 03 the zero set in p2(K) of {pi,...,p4}. Then, the al- 
gorithm decomposes P^(K) \ ® as 

P2(K) \ 23 = uLi©fc 

such that if (to : 5o : Do) ^ ©fe then P(to,5o, Dq) is a point of S of multiphcity k. 

As already remarked in Section [2], we assume that none of the polynomials pj is 
zero or, more generally, that there do not exist pj,pj and A G IK such that pj = Xpj. 
Note that this excluded situation corresponds to a plane, and hence ©i = P^(K). 

In addition, we use the notation introduced in Section [21 namely the affine surfaces 
iSa;-, the affine rational parametrizations Pa;-(t, s), and the polynomials Pk,<li,j- More- 
over, we also use the notation Vxi{i,v), Vxi{s,v) (see Section E]). In this situation, the 
algorithm is as follows. 

Algorithm 

[Preparatory Steps] 
[Step 0.] If any of the projective curves defined by the non-constant polynomials 
in {pi,P2,P3,5'4,i, '?4,2, ^4,3} of the parametrization Vx4^ passes through (0 : 1 : 0) we 
perform a suitable polynomial linear change in the parameters {t, s}. 

[Step 1.] Apply Method 1 to compute 722 := degiV^J (see Section [3]). 

[Step 2.] Apply Method 4 to determine an affine point, say Aq, out of the affine surface 
Sx^ (take Aq with non-zero components such that if the algorithm, in subsequent steps, 
requires a point in K^ \ Sx^ with z 7^ 4 no further computation would be needed (see 
Remark [3l3l) and apply Method 2 to compute ni := $2,3(^0) (see Section [3]). 
[Step 3.] Let Aj := {(to,so) \lcm{qi^i, qi^2, qi,3)(to, sq) = 0}, z = 1, . . . ,4. 
[Step 4.] Let ^^ = flti ^i (see Lemma [22D find 03. 

[P-affine singularities (First level: part I)] 

[Step 5.] Compute (see Section [1 A* = Vx,{t*,s*); ^2,3{A*) = (xf ,xf ); Q^''"^* = 
Numer(xf (t,s,r,s*) - xf {h,h2,t*, s*)); as well as K = gcdiQ^^''^\Q^''''-^*). 
[Step 6.] If deg|4,|(A') > then return ©1 = p2(K) {S is a plane). 
[Step 7.] Computation of Aq 
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[Step 7.1.] Compute the leading coefficient Tj of 0^^'*' w.r.t. t, 

[Step 7.2.] Compute r := Resj(0i''''' ,52"^* )• 

[Step 7.3.] Zx is the zero set of at least one non-zero coefficient, w.r.t. {hi, /i2}, of 

the homogeneous form of maximum degree of K. 

[Step 7.4.] Aq is the zero set of all coefficients of Ti, T2 w.r.t. {hi, /i2, s} union the 

zero set of all coefficients of r w.r.t. {hi, /12, s} union Zk- 
[Step 8.] Computation of A^ 

[Step 8.1.] Homogenize (w.r.t. {t, s}) the numerators and denominators of xf*: ^-^d 

substitute them in (0 : 1 : 0). Take Zoo as the union of the zero sets in K^ of these 

polynomials. 

[Step 8.2.] Compute Hi = Numer(x^ '''*' — Xj) as well as the leading coefficient Mi 

of Hi w.r.t. s. 

[Step 8.3.] A^ is the zero set of all coefficients of Mi w.r.t. {t,Xi,X2} union the 

zero set of all coefficients of M2 w.r.t. {t, Xi,X2} union Zoo- 
[Step 9.] Computation of Aj 

[Step 9.1.] Compute R := ReSs{Hi, H2) and its leading coefficient W w.r.t. t. 

[Step 9.2.] A2 is the zero set of all coefficients of W w.r.t. {Xi,X2}. 
[Step 10.] Computation of A3 and A4 

[Step 10.1.] Compute the set {Uij \ {i,j) G J} of all coefficients of i? w.r.t. {Xi,X2}. 
Step 10.2.] Compute a = gcd({ajj | {i,j) G J}) and ajj = ^. 

Step 10.3.] Determine the leading coefficient N oi a w.r.t. t. 

Step 10.4.] If Card(J) = 2 (say J = {(zq, Jo), (^1, Ji)}) 

[Step 10.4.1.] Compute the leading coefficient Lj of ca^jj w.r.t. t (j = 0, 1) and 

[Step 10.4.2.] A3 is the zero set of {Lq, Li} union the zero set of S. 
[Step 10.4.3.] A^ is the zero set of N. 
[Step 10.5.] If Card(J) > 2 (say J = {iik,Jk)}k=i,...,e) 

[Step 10.5.1.] Compute the leading coefficient L of aj^j^ w.r.t. t and S = 



e 



ReS((ai^j^, ai2j2 + X]fc=3 W^fc-2«i 



k iJk ) 



[Step 10.5.2.] A3 is the zero set of all coefficients of S w.r.t. \Wi, . . . , Wi-2}- 
[Step 10.5.3.] A4 is the zero set of L union the zero set of A^. 
[Step 11.] Set A* = UtoA*, and include i{(K^ \ A4) \ A*) in ©i. 

[P-affine singularities (First level: part II)] 

[Step 12.] Decompose A* into irreducible components. 

[Step 13.] For each point A at a zero- dimensional component of A*, HA ^ A4 then 

apply Method 5 to compute a = Tault{A,Sx^), and include }{A) in ©„. 

[Step 14.] For each 1-dimensional irreducible component C of A*, compute its genus. 

[Step 15.] If C is rational proceed as in Section Hj (Level 1, Part II, rational case). This 

will generate an open subset C* of C where the multiplicity is invariant and that would 

be included, via j, in the corresponding &k- For the finitely many points in the C \ C*, 

proceed as in Step 13. 

[Step 16.] If C is not rational proceed as in Section [4] (Level 1, Part II, positive genus 
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case). This will generate an open subset C* of C where the multiplicity is invariant and 
that would be included, via j, in the corresponding (3^. For the finitely many points 
in the C\C*, proceed as in Step 13. 

[P-affine singularities (Second, Third and Fourth Level)] 
[Step 17.] If A4 \ ^a = (see Step 4) go to Step 20 else proceed as follows 



[Step 17.1 
[Step 17.2. 
[Step 17.3. 

[Step 18.] If 
[Step 18.1. 
[Step 18.2. 
[Step 18.3. 

[Step 19.] If 
[Step 19.1 
[Step 19.2 
[Step 19.3 



If A4 \ A3 = go to Step 18. 
Compute the irreducible decomposition of A4 \ A3. 
Proceed as in Steps 13, 14, 15, using "Pj-g instead of "P^^^. 
(A4 n A3) \ QSq = go to Step 20 else proceed as follows 
If (A4 n A3) \ A2 = go to Step 19. 

Compute the irreducible decomposition of (A4 fl A3) \ A2. 
Proceed as in Steps 13, 14, 15, using Vx2 instead of Px4- 



(A4 n A3 n A2) \ 5Sa = go to Step 20 else proceed as follows 
If (A4 n A3 n A2) \ Ai = go to Step 20. 

Compute the irreducible decomposition of (A4 fl A3 fl A2) \ Ai. 
Proceed as in Steps 13, 14, 15, using Vx^ instead of Vx4^- 

[P singularities at infinity] 
[Step 20.] If (0 : 1 : 0) ^ 03 (i.e., not all pi(0,l,0) vanish) apply Method 5 to com- 
pute a := m.u\t{Vx4^{0,0),SxJ and include (0 : 1 : 0) in ©q,; we are assuming that 
p4(0, 1, 0) 7^ 0, otherwise take other component and proceed accordingly. 
[Step 21.] Check whether p4(l,A,0) does not vanish. If it does vanish, find pj not 
vanishing at (1 : A : 0) and proceed accordingly. 

[Step 22.] For each (1 : Aq : 0) such that p4(l, Aq, 0) = if (1 : Aq : 0) ^ ®: find pj such 
that pj(l, Ao,0) 7^ 0, compute a := m.Vi\t{Vxj{Xo,0),Sxj), and include (0 : 1 : 0) in &a. 
[Step 23.] Proceed as in Step 15 using 'Pj.^(s,-0), instead of P^^^, and the curve (A,0). 

Example 6.1. We consider the parametrization 

P(t,5,D) = (s^ iS^ + t^ + D^ : (t+2s)D : (s + t) t)) 

of the surface S. One can easily check that the parametrization satisfies all hypotheses 
in Section [2] In addition 



^x, 



p 



s^ + r + 1 t + 2s s + t 



s"^ ' s^ 



t + 2s 



V 



s^ + r + 1 s + t 



t + 2s t + 2s t + 2s 



s + t 



s^ + t^ + V s^ + t^ + V s^ + t^ + 1 



V 

I X4, 



-2 I j-2 



t^ + 1 t + 2s 



S + t' s + t 's + t 



Note that P^^ satisfies the hypotheses in Step 0. In Step 1 one gets ?t,2 := 1, and in 
Step 2 we get Aq := (1, 1, 1) and rii := 4. In Step 3 we get that 

• Ai is the line s = 0, 
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• A2 is the complex circle s + t + 1 = 0, 

• A3 is the line t + 2s = 0, and 

• A4 is the line s + t = 0. 

Therefore, in Step 4 we get Q3a = 25 = 0. In Step 5 we get that K = t* + s*. We start 
the computation of Aq. In Step 7.1. we get 

Ti = (S* + t*) {-h2^S* - h^h* + 3*^2 + S*\i) , T2 = (s* + t*) h2 s* {-s* + /la) . 

In Step 7.2. 

r = -[s* + t*f{s - s*){-s + h2){-hls* - hit* + s*^h2 + s*^h^){-hi^s*'^s + /ii2s*2s2 + 
s*^hi^sh2 - 2 /i2 r s^s* /zi + s*^h2 - 2 /z^s*^ + r^sh^s* + his* + sH*^hl - st*^hl). 
In Step 7.3. Zk is the line s* + t* = 0. Finally, in Step 7.4. we conclude that 

A* = {(t*,s*) GK^ir + s* = 0.} 

In Step 8. we compute A^. In Step 8.1. we get that Z^ is union of the lines t* + s* = 
and t* = 0. In Step 8.2. we get M^ = s* + 1* - Xi(s* + 1*), M2 = -X2(s* + 1*). So, in 
Step 8.3. we conclude that 

A* = a; U {(t*, s*) G K^ 1 1* = 0}. 

For the computation of Ag, in Step 9.1. we get W = X.l{s* + t*)^. Therefore, 

A* — A* 

For computing A3 and A4, in Step 10.1., we get that R has 6 non-zero coefficients 
w.r.t. {Xi,X2}. Moreover, a = {s* + t*f{t - t*) (see Step 10.2.) and N = {s* + t*f 
(see Step 10.3.). Since Card(J) = 6 we go through Step 10.5. Then (see Step 10.5.1.), 
ai^j^ = t{s*^ + t*2) -t*, L = s*^ + t*^, and 
S = -t* s*^ - 2 t*^s*^ - t*^s*^ + 2Wi s*H*^ - Wi t*^s*^ + 2 s*'^W2 1*^ + 2 s*^W2 1*^ - 

2 s*^W2 1*^ -2W4 s*H*^ + W4 t*h*^ + W4 s*H* + 2 t*^Wi s*^ + 2 t*^W4 s*^ + 2 s*^Wi t* - 

3 s*^Wit*^ -3 s*^Wit*^ + 5 s*^W2t*'^ + 4:S*^W2t*'^ -2 s*^W2t* + s*'^W3t* + 2 s*^W3t*^ - 

4 s*^W4 1*^ + 3 s*H*^Wi + 5*^1^4 1* + 3 s*H*^Wi - t*^Wi s*2 + 2 t*^s*^W2 + t*^s*^W^ - 
2 t*^Wi s*^ + t*'^W4 s*^ + s*^Wi - s*^Wi + s*^W2 - 2 s*'^W4 

In Step 10.5.2. and Step 10.5.3. we get 

A* = {{t*, s*) eK^\s* = 0}, A* = {(r, s*) e K^ \ s*^ + t*^ = 0} U A*. 
Therefore in Step 11. we get 

A* = {(t*, s*) G K^ I s* = 0} U {(t*, s*) G K^ 1 1* = 0} 

U{(r, s*) G K^ I s* + r = 0} U {(t*, s*) eK^\s* ± it* = 0} 

Now Step 12. is already executed. Step 13. is not needed, and in Step 14. we get that 
all components are rational; indeed lines. Then, for each of the lines we execute Step 
15. 
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Let C be the line s* = 0. We consider the normal proper parametrization 
Q(A) = (A, 0), and H(A) = (o, i^, i^J . 

C n Aq = {(0,0)} C A4. So, we deal generically with 'H(A), knowing that 
g7'^4,H(A) ^ i_ rj.^^^^ jg^ ^g g^ ^^^j^ ^^ g^gp g taking ^* as 7/(A). We get that 

K = 1. We know that A;^ = 0. In Step 8. we get that the new Z^ = {0} C K, 
that Ml = A(l - Xi),M2 = -X2. Therefore, the new Al = {0} C K. In Step 
9. R = XjXt^ + XlX - Xp - X^A^t - X2 - X2A2 + A - Xi and VT = X^A. 
So A2 = A^. In Step 10.1., we get that R has 4 non-zero coefficients w.r.t. 
{Xi,X2}. Moreover, a = 1 (see Step 10.2.) and A^ = 1 (see Step 10.3.). Since 
Card(J) = 4 we go through Step 10.5. Then (see Step 10.5.1.), flj^jj = A, L = A, 
and 5 = A^. Thus, A3 = A4 = A2 = AJ. Summarizing, for all A 7^ we 
deg($2,3(^(A)) = degi(PrimPart{Xi,x2}(Res,(iJi,/72))) = 2. Therefore, rm = 2 
and then 

mult((A : : 1)) = muh(?/(A),5^J = ^^ ~ "^^ = 2, with A ^ 0. 

n2 

Let C be the line t* = 0. We consider the normal proper parametrization 

Q(A) = (0, A), and ?/(A) = (\, ^^, 2) . 

C n Aq = {(0,0)} C A4. So, we deal generically with 'H(A), knowing that 
gP^4,-H(A) ^ i_ rjj^g^^ jg^ ^g g^ ^^^^ ^^ g^gp g taking A* as 'H(A). We get that 

K = 1. We know that Aq = 0. In Step 8. we get that Zoo = K. So, we 
perform a suitable linear change of parameters in Vx^ to avoid that, namely 
we replace (during the analysis of this curve) Vx4^ by Vx4^{s — t,t + s). Then, 
we get that Z^o = {0}. Then, repeating the computation we get that A3 = 
A* = a; = A* = {0}. Summarizing, for all A 7^ we deg($2,3('H(A)) = 
degj(PrimPart{Xi,X2}(R'eSs(i/i, i/2))) = 3. Therefore, mi = 3 and then 

muh((0 : A : 1)) = mult(?{(A),5^J = I^IZJIH = 1^ with A ^ 0. 

n2 

The next curve is precisely A4. So, we postpone its analysis to further levels. 

Let C be the lines t* ± i s* = 0; we treat both curves simultaneously. We consider 
the normal proper parametrization 

Q(A) = (±^A,A), and?/(A) ' ^ ^ ±«A + 2A 



A± tX X±tX X±tX 
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C n Aq = {(0,0)} C A4. Repeating the computation we get that A3 = 
Al = A*2 = Al = {0}. Summarizing, for all A 7^ we deg(<l>2,3('H(A)) = 
deg^{PnmPait^Xi,X2}{^^^s{Hi,H2))) = 3. Therefore, mi = 3 and then 

mult((±zA : A : 1)) = mult(H(A), 5^J = ^^ ~ ^^ = 1, with A 7^ 0. 

We go to Step 17. A4 fl A3 = {(0,0)}. So, we work generically with A4 and Vx^- 
So we consider Q(A) = (—A, A) and 'H(A) = Vx^iQiX))- Proceeding as above, we 
get A*2 = Al = A*2 = Al = {0}. Summarizing, for all A 7^ we deg($2,3('H(A)) = 
degj(PrimPart{Xi,X2}(R'6Ss(iJi, 7/2))) = 3. Therefore, nii = 3 and then 

mult((-A : A : 1)) = mult('H(A), 5^J = ^^ ~ "^^ = 1, with \ ^ 0. 

n2 

In Step 18. since (A4 fl A3) \ A2 = {(0, 0)}, we compute the multiplicity of 7^x2(0, 0) = 
(0,0,0) using P^^a- We get 

mult((0 : : 1)) = mult((0, 0, 0),5^J = 3. 

Since A4 fl A3 fl A2 = we skip Step 19. and we pass to Step 20. 
In Step 20 we first observe that ^ = 0. Moreover, since p4(0, 1,0) = but 
pi(0, 1,0) 7^ we compute mult('Pi.^(0, 0),iSi.J by applying Method 5. One gets that 
mi:=deg(<l>2,3(Pxi(0,0)) = 2. So, 

mult((0 : 1 : 0)) = mult(P^JO, 0),5^J = ^^ ~ "^^ = 2. 

n2 

In Step 21. we observe that p4(l. A, 0) = but pi(l. A, 0) = A^. In Step 22. we need to 
analyze (1:0:0). We do it using Vx2{s,v) to get 

mult((l : : 0)) = mult((0, 0, 0), 5^^ = 3. 

In Step 23. working with Vxi{s,v) we conclude that 

mult((l : A : 0)) = mult ( [ ^^, 0, 0), SxA j = 2 for A ^ {0, 1, -1} 

So it only remains to analyze (1 : 1 : 0), (1 : — 1 : 0). We apply Method 5 with Vx^ to 
get 

mult((l : 1 : 0)) = l,mult((l : -1 : 0)) = 1. 

In Fig. 1, we summarize the conclusion. 
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Line at infinity (t:s:0) | 



(1:a:0) (1:0:0) 
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Affine plane (t:s:1) | 
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(0:s:1) 



,(0:0:1) 
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(t:0:1) 



Figure 1: Decomposition of the parameter space P^(K), where the multiphcities 
mult((t : s : v)) are represented instead oi {t : s : v) 
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